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In a companion paper we have introduced a model of scalar field dark energy, Cuscuton, which 
can be realized as the incompressible (or infinite speed of sound) limit of a k-essence fiuid. In this 
paper, we study how Cuscuton modifies the constraint sector of Einstein gravity. In particular, we 
study Cuscuton cosmology and show that even though Cuscuton can have an arbitrary equation 
of state, or time dependence, and is thus inhomogeneous, its perturbations do not introduce any 
additional dynamical degree of freedom and only satisfy a constraint equation, amounting to an 
effective modification of gravity on large scales. Therefore, Cuscuton can be considered to be a 
minimal theory of evolving dark energy, or a minimal modification of a cosmological constant, 
as it has no internal dynamics. Moreover, this is the only modification of Einstein gravity to 
our knowledge, that does not introduce any additional degrees freedom (and is not conformally 
equivalent to the Einstein gravity). We then study two simple Cuscuton models, with quadratic and 
exponential potentials. The quadratic model has the exact same expansion history as ACDM, and 
yet contains an early dark energy component with constant energy fraction, which is constrained to 
^ 2%, mainly from WMAP Cosmic Microwave Background (CMB) and SDSS Lyman-a forest 
observations. The exponential model has the same expansion history as the DGP self- accelerating 
braneworld model, but generates a much smaller integrated Sachs- Wolfe (ISW) effect, and is thus 
consistent with the CMB observations. Finally, we show that the evolution is local on super-horizon 
scales, implying that there is no gross violation of causality, despite Cuscuton's infinite speed of 
sound. 



I. INTRODUCTION 

The nature of the current acceleration of cosmic expan- 
sion is among the most outstanding puzzles in theoret- 
ical physics. Various cosmological observations such as 
the dimming of distant supernovae la 0, @| > anisotropies 
in the Cosmic Microwave Background (CMB) M, and 
the large scale structure of the Universe (e.g. 0|) can 
be most easily explained by having an exotic dark energy 
component with negative, nearly constant and uniform, 
pressure which dominates the energy density of the Uni- 
verse. 

While the simplest model of vanilla dark energy, i.e. 
a cosmological constant, remains consistent with all the 
present observations many models of non-minimal 
dark energy have been developed in anticipation of any 
future failure of cosmological constant in explaining the 
data. However, all the models that predict an observable 
dark energy density evolution suffer from an extreme fine 
tuning problem, requiring incredibly light masse scales 
(^ 10^'^^eV) that are hard to protect from quantum cor- 
rections. 

In a companion paper [l[ , we developed a new model of 
field theoretical dark energy, Cuscuton, that while gener- 
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ally non-uniform, lacks any dynamical degree of freedom. 
This protects the theory from quantum corrections at 
low energies, and thus makes it a perfect candidate for 
an evolving dark energy in the current era. The name 
Cuscuton (pronounced kas-kii-tan), is derived from the 
Latin name for the parasitic plant of dodder, Cuscuta. 
Classically, it is a new kind of constraint system, allow- 
ing a novel class of constrained dynamics. 

The Cuscuton action for the scalar filed (p can be writ- 
ten as 

S = J d^xV^[^,^^\g^-^d^^d.^\ - U(^)], (1) 

where fi is an (arbitrarily defined) energy scale. One 
can show that Cuscuton is an incompressible k-essence 
fluid @, 0]j he. it has an infinite speed of sound, raising 
questions about the causality of the theory. However, in 
we show that there is no breakdown of causality, as 
the perturbations lack (symplectic) dynamics, and thus 
do not carry information. 

Given that Cuscuton acts as a constraint system, it is 
interesting to ask how a minimal coupling of Cuscuton 
to gravity will modify interaction of gravity with matter. 
In particular, as is well known, the scalar perturbations 
about the FRW metric represent the general relativistic 
spacetime constraints for matter coupled to gravity. A 
non-trivial choice of Cuscuton potential will modify the 
constraint equations, leading to observable consequences 
in cosmology. We aim to characterize these observables 
associated with the scalar perturbations. 

We derive general scalar perturbation equations in the 
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presence of Cuscuton and show analytically and numer- 
ically how the constraints are modified. We show that 
Cuscuton models can have an expansion history identi- 
cal to that of ACDM but with different CMB and mat- 
ter power spectra due to gravitational potential evolu- 
tion mimicking that of tracker models of quintessence. 
We also find that Cuscuton can exactly replicate the ex- 
pansion history of the DGP self-accelerating cosmology, 
while predicting similar small angle CMB and matter 
power spectra, at a few percent level. The main differ- 
ence between the DGP self-accelerating cosmology and 
Cuscuton will show up due to the anisotropic-stress in- 
duced boost in the Integrated Sachs- Wolfe (ISW) effect 
for the DGP model. Finally, we briefly consider the issue 
of causality in the scalar sector and find that a Yukawa- 
like exponential fall-off of Cuscuton effects beyond the 

scale protects the gravitational phenomenology from 
gross anomalies in the long wavelength limit. We show 
that the usual Bardeen parameter remains approximately 
constant in the long wavlength limit. 

The paper is organized as follows: In Section |TT1 we 
compute the background evolution, followed by linear 
perturbations in Sections IIIII We then study how the 
matter power spectrum and the ISW effect are affected by 
Cuscuton'm Section llVl and discuss current observational 
constraints. Before presenting our conclusions in Section 
IVIl we also discuss causality of linear Cuscuton pertur- 
bations in Section |Vl Throughout this paper, we use the 
reduced Planck's constant Mp = {SttGn)'^^'^ ~ 2 x 10^^ 
GeV. 



II. BACKGROUND EVOLUTION 

Starting with the Cuscuton action (Eq. [T]) in a flat 
Friedmann-Robertson- Walker (FRW) metric 



ds^ = dt^ - a{tfdx'dx' 



(2) 



and a homogenous field configuration, i.e. Lp ~ '^{t), the 
action takes the form 



a^dt[^l^\^\-V{v)\ . 



(3) 



Varying the action with respect to yields the back- 
ground field equation 



{?,^l^H) sgn(^) + V'{^) - 0. 



(4) 



One implication of Eq.Q is that V{(p) always de- 
creases (increases) with time in an expanding (collapsing) 
universe. Eq.Q, in combination with the Friedmann (or 



Gq Einstein) equation in a flat universe 



Ptot 

3M2 



yields 



(5) 



(6) 




FIG. 1: Different possible Cuscuton potentials in the V — V' 
phase space. The shaded region is excluded in a flat FRW 
cosmology, as it requires negative matter density, p,n (Eq. 
|6|. The dark solid, dashed and dotted lines are quadratic 
potentials with zero, positive and negative bare cosmological 
constants. The red (dark) arrows show the time evolution of 
the field in a flat universe, and the points show the asymptotic 
state of the field, if present. The blue (light) solid line shows 
a constant pm > contour. As we show in the text, in addi- 
tion to avoiding the excluded region, the potentials should be 
shallower than the constant pm contour in the V — V' plane. 



where pm is the background density of ordinary matter 
in the Universe. Notice that the Cuscuton kinetic term 
does not contribute to its energy density. 

In Fig.ll]), we show different possibilities for the Cus- 
cuton potential in the phase space. Eq. ([5]) implies 
that for a positive matter density {pm > 0; weak energy 
condition), assuming a flat universe, we should have 



n.,<(0- 



(7) 



The excluded region is shown by the blue shaded area 
in Fig. Il]), while the blue (light) line shows a constant 
Pm > contour. 

The continuity equation for matter density pm — 
—3H(pm+Pm), in combination with Eq.Q and the time 
derivative of Eq.®, assuming the null energy condition 
Pm + Pm > 0, gives: 



2M 



(8) 



which puts an additional constraint on the potential: 



2 dV{ip) 



(9) 



Therefore, all the allowed potentials need to be shallower 
than the constant density contours (blue/light line) in 
Fig.©. 
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Dark lines in Fig.(IT]) show different quadratic poten- 
tials of the form 



(10) 



Interestingly - as we will show in Sect ion HVl - a quadratic 
potential leads to a tracking behavior. In this case, the 
quadratic term of the potential always maintains a con- 
stant fraction of the total density of the Universe. The 
dark solid line shows a quadratic potential with no bare 
cosmological constant (Vq = 0). For a positive bare cos- 
mological constant (Vb > 0; dashed line), it is interesting 
to notice that the potential never reaches its minimum, 
and is stalled at the boundary of the pm < region 
(shaded area). Therefore, the effective value of the cos- 
mological constant is larger than its bare value. Finally, 
a potential with a negative bare cosmological constant 
(Vo < 0; dotted hne) results in a bounce, as V , and thus 
the expansion rate go to zero, and changes sign subse- 
quently. 

Another interesting example is an exponential poten- 
tial of the form 



V{ip) = Vo exp 



(11) 



Direct substitution in Eqs. yields 



H = 




(12) 



which is exactly the same as the background dynamics 
in a (flat) Dvali-Gabadadze-Poratti (DGP)[i] 5D self- 
accelerating brane- world model Q . However, this is only 
a coincidence, as the detailed dynamics of metric pertur- 
bations cannot be identical. For example, the anisotropic 
stress always vanishes in Cuscuton models, which fol- 
lows from the use of 3-l-lD Einstein equations for lin- 
ear perturbations of a scalar field. On the contrary, the 
anisotropic stress is generically non-vanishing (see e.g. 
[lOj and references therein) in DGP models. Neverthe- 
less, this coincidence can be used to examine the ob- 
servable differences between dark energy and modified 
gravity models with exact same background dynamics. 

Before concluding this section, it is important to em- 
phasize that Cuscuton is classically a theory of modi- 
fied gravity, rather than Einstein gravity with additional 
scalar field degree of freedom. This can even be seen 
at the level of background equations, as combining Eqs. 
(|3]in|), (for ip < 0) we can write the modified Friedmann 
equation as 



+ y [y'-i(3M'i7)]}, 



(13) 



where 1/'"^ is the inverse function of V'{ip). We see that, 
in general, is no longer linearly dependent on the en- 
ergy density pm and the exact nonlinearity is controlled 



by the choice of function V{lp). Moreover, unlike in ordi- 
nary Einstein-Hilbert action coupled to a homogeneous 
scalar degree of freedom, the modified Friedmann equa- 
tion is fixed once p,n is fixed, i.e. one does not need to 
specify initial/boundary conditions for Lp. 

To obtain intuition for how the modification works, 
let us again consider the quadratic potential. If V{(p) — 
^m?ip^ , it is simple to check that Eq. (|13p is identical to 
the ordinary Einstein-Hilbert Friedmann equation with 
a renormalizcd Planck's constant 



Mi 



Mi 



2inn? 



(14) 



This is a manifestation of Cuscuton modification of grav- 
ity. 

This renormalization of the Planck mass is reminiscent 
of a similar effect for Lorentz- Violating vector fields [ll[ , 
as well as the exponential quintessence model (see Sec. 
lIVBp . although, in contrast to Cuscuton, both models 
do introduce an additional dynamical degree of freedom. 



III. LINEAR PERTURBATIONS IN CUSCUTON 
COSMOLOGY 

A. Linearized Field Equation 

Varying the Cuscuton action with respect to in a 
general curved space-time yields 



X 



V^V> + /i~2\/XU'(v3) 0, (15) 



where X = d'^(pd\(p, and denotes covariant deriva- 
tive. Assuming a linearly perturbed FRW metric in the 
longitudinal gauge 



ds^ = (1 + 2<i>)dt^ - a{tf{l - 2^)dx^dx' 



(16) 



we can evaluate the field equation at the linear order in 
field/metric perturbations: 

3ip{^ + H<P) + a-^V^Sip - p-^\ip\V"{ip)Sip = 0. (17) 

Here, V'^ is the spatial Laplacian with respect to comov- 
ing coordinates. It is interesting to note that the linear 
perturbation equations do not include any second order 
time derivative. 

Taking the time-derivative of Eq. we find 



p^ 



-3H. 



(18) 



This yields the following form for Eq. ([TT]) in the Fourier 
space: 



3(^($ -I- H$) 



(19) 
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explicitly showing that, as pointed out in Cuscu- 
ton perturbations simply follow metric perturbations in 
a non-local way, and do not introduce any additional dy- 
namical degree of freedom. In other words, similar to the 
homogenous equation (j4]) , the field equation only amounts 
to a constraint condition that relates metric and field per- 
turbations. 



B. Linear Einstein Equations: Modified Gravity 



The Einstein equations for scalar metric perturbations 
in the presence of Cuscuton as well as ordinary (dust) 
matter inhomogeneities are simply written as: 



$ + H$ = (2M2)-1(^2^(^ ^ 



(20) 
(21) 
(22) 



where A is the potential of the matter peculiar velocity ~ a^^dX/dx^. 

Transforming to the Fourier space, it is interesting to note that dip can be eliminated by combining the field 
equations (Eqs. [4l and fT9|l and the Goo equation (Eq. [20]) to yield a modified law of gravity 



3H 



9H{2H + SH^nm) 



($ + i/$) + {2M^)-^6prn = 0, 



(23) 



where il^ ~ Pm/ (■^MpH^), is the matter density in units of the critical density of the Universe. 



Thus we notice that Sip completely drops out of the 
linear gravity (or Poisson) equation, although it modifies 
the equation in a non-local way. This is another mani- 
festation of Cuscuton being a theory of modified gravity, 
even though it is a particular limit of k-essence. 

It is important to note that this modification (or 
screening of gravity) does not affect Newtonian gravity 
on small sub-horizon scales, i.e. as long as ^ ^ H , H. 

We should point out that the only other modification 
of Einstein gravity that does not introduce a new de- 
gree of freedom (to the best of our knowledge), known 
as Modified-Source Gravity [l3|, is conformally equiva- 
lent to the Einstein gravity. Moreover, Modified-Source 



Gravity can be constructed by a non-linear local modifi- 
cation of the matter Lagrangian, where, in terms of the 
modified Lagrangian, (and in contrast to Cuscuton) the 
gravity reduces to Einstein gravity. 



C. Linear Einstein Equations: Modified Dynamics 

Let us consider the modified dynamics of the gravi- 
tational potential, in the presence of Cuscuton and 
pressureless dark matter. Plugging Eq. ([19]) into Eq. (|22p . 
after straightforward manipulations, we arrive at 



(1 + C2)$ + (4i? + Ci + C2H + C3)$ + {3H'^ + H- -VL„,H'^ + CiH + C^H + C^H)^ = 0, 



^1 = — p — Z7Z — )'^2 = — — ,03 = — 



2H{!^Y+3H {2H + 3H^nm) 



3H 



r 



(24) 
(25) 



which is the desired differential equation for $. Note 
that in the limit of pure matter dominated case without 
Cuscuton modification of gravity, we have 



implying that the solution to Eq. ([24]) asymptotically 
approaches $ — >constant. Also, note that any nontrivial 
scale dependence introduced by Cuscuton is character- 
ized by the scale H, and not H. 



1, and 2H + 3H^,C^ 0, 



(26) 
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Consider the long wavelength limit {k/a)"^ <C H- First, 
let us ask whether a constant $ is a solution to linearized 
Einstein equations with just pressureless dust field de- 
grees of freedom. The Cuscuton modification that could 
prevent $ from being a constant is the coefficient of $ in 
Eq. 



(27) 



3H^ + H- -n„^H^ + CiH + C2H + C3H 



HH 
2IP 



IV. OBSERVATIONAL SIGNATURES OF 
CUSCUTON COSMOLOGIES 



In this section, we study the observational signatures 
of Cuscuton cosmology. First we consider the observables 
analytically, in a perturbative setting. Afterwards, we fo- 
cus on different observational signatures and constraints 
for quadratic and exponential Cuscuton potentials. 
A. Analytic Treatment of General Cuscuton 
Potentials 



which generically does not vanish unless the scale factor 
is of the form 



a = a,:[l + ^(t- 



ti 



]Hf/\H. 



(28) 



corresponding to power-law expansion, or a constant ef- 
fective equation of state. Therefore, having a constant 
effective equation of state will allow a constant <i> solu- 
tion on long wavelengths. 

Next, whether or not $ is damped (due to the friction 
term in Eq. [M)) depends on the sign of 



4:H + Ci+ C2H + C3 
I + C2 



H 



dt 



InH- 



(29) 



Hence, we see that whether the potential decays or not 
depends on how fast H horizon grows/decays with time. 
In the case a constant effective equation of state (or Eq. 
\28\i . the potential will always decay until it asymptoti- 
cally reaches a constant. As we will argue in Section Ivl 
the general behavior of <I> on superhorizon scales can be 
easily understood from the conservation of the Bardeen 
parameter. 

In the opposite limit of (fc/a)^ ^ H, the coefficient in 
Eq. that prevents $ from being constant becomes 



3i/^ 1 - 



H 



(30) 



which again does not necessarily vanish with the grav- 
ity modified by a nontrivial V(ip), even though the form 
of the terms looks identical to that of Einstein gravity 
sourced by pressureless dust. That is because the back- 
ground Einstein equations are modified by the presence 
of Cuscuton. Note that, as far as the damping coefficient 
is concerned, since Ci — > in the short wavelength limit, 
the $ term always acts as a damping term. 



In addition to the matter power spectrum, the decay- 
ing gravitational potential caused by the Cuscuton mod- 
ification will manifest itself in the anisotropy of the Cos- 
mic Microwave Background (CMB). The induced CMB 
anisotropy due to the Fourier mode $k is til] : 



rvo 

e,,k = / d?7<?(77)(eo + $k)ji[fc(?7o -??)] + 
— J dri Vb{k)g{r])—ji[k{f]o - f])] + 
2 d^^e-^^'^^mm^v)], (31) 



where rj is the conformal time, Qq is the monopole tem- 
perature fluctuation, 3(77) is the visibility function, Vb is 
the baryon velocity field, r is the optical depth, and 770 is 
the present day conformal time. The first term is known 
as the Sachs- Wolfe effect, while the second term is the 
dipole contribution, and the last term is the Integrated 
Sachs- Wolfe (ISW) contribution, which depends on the 
time variation of the gravitational potential. All the in- 
tegrals are taken over the light cone. 

We can compute the ISW contribution to CMB 
anisotropics perturbatively. To first order in the Cus- 
cuton potential, V, one can write d^/drj, which appears 
in the ISW term of Eq. as 



(l).k 



drj 



u 



t(v) /f' 



8/3 



Skit'), (32) 



where t{r]) '-^ a^rf /27t~, S{t) is given by 



^k(i) 



2,tV 



4 id + ^ 



V V 
I 



2t 



(0),k 



3M2 



(33) 



where '&(o),k is the zeroth order solution in V, to Eq. (|24p . which is a constant in a flat matter dominated uni- 
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verse. The key feature manifested in Eq. (|33|) is the mild- 
ness of the scale dependence. This is simply a result of 
the fact that Ci ^ in the limit of fc oo, while Ci 
does no vanish in the limit of /c — > 0. Although more 
explicit expressions for O; and Ci may be derived, gen- 
eral results deviating strongly away from cosmological 
constant-induced ISW effect are complicated and unillu- 
minating. Hence, in the subsequent subsections, we will 
numerically examine two interesting explicit potentials, 
and compare them against current cosmological observa- 
tions. 



B. Quadratic (Tracking) Potential 



Using potentials of the form V{ip) 
one finds solutions similar to that of a tracking dark en- 
ergy component plus a cosmological constant. 
While Vo simply contributes towards the cosmological 
constant, the quadratic term ^rn^Lp^ maintains a con- 
stant fraction of the total energy density of the Universe, 
much like models of early dark energy [H, • To 
see this, take the square of Eq. ji]) together with the 
Friedmann Equation (O to get 



AMI 
M2 



1 ? 9 
Plot 



2M2m2 



const., 



(34) 



where AM^ denotes the equivalent change in the large 
scale Planck mass (Eq. [H)) . Therefore, assuming Vq ^ 0, 
the expansion history H(z) is exactly equivalent to that 
of a ACDM cosmology, as the extra quadratic compo- 
nent simply follows the rest of the energy content of the 
Universe. In fact, the expansion history H{z) becomes 
independent of Qq , provided one scales the other compo- 
nents of the Universe appropriately, i.e. simultaneously 
transforming fi^ — > {l — il.Q)^i, where i labels A, photons, 
baryons and dark matter while changing the number of 
relativistic neutrino species N^, = 3.04(1 — 2. 47r2Q) leaves 
H{z) invariant. This leads to two corollaries: First, there 
will be a degeneracy between and Jig (at least, as long 
as we only consider the expansion history). Secondly, 
quadratic Cuscuton can only be constrained by studying 
sub-horizon perturbations, provided one does not a priori 
know the number of neutrino species. 

Current constraints on light element abundances, 
which are predicted from Big Bang Nucleosynthesis, al- 
ready put significant constraints on the number of rela- 
tivistic neutrinos during the radiation era, — 3.1±0.7 
(2Q1] . which yields Q.q ^ 10%. However, as we see be- 
low, observational constraints on cosmological inhomo- 
geneities can put much tighter constraints on ilg. 

Since Cuscuton perturbations have an infinite speed 
of sound, they do not cluster on sub-horizon scales, and 
thus perturbations start to decay as they enter the hori- 
zon. This behavior is reminiscent of scalar field dark en- 
ergy (or quintessence) for which, once inside horizon, the 
perturbations also free-stream with the speed of light. 



One can work out this sub-horizon decay analytically 
during matter domination where a oc ^2/3. Then, Ignor- 
ing 6ip in Eq. (gj) (or using Eq. [^Iwith d ^ 0), in the 
k ^ aH regime, we end up with: 



$ + 4ff<i> + = 0, 

which can be easily solved ii H — 



$ cx a";Q! = - 
4 



-1 ± 




(35) 



(36) 



This result is identical to the growth suppression for a 
quintessence field with an exponential potential, which 
has a similar tracking behavior [45[ (e.g. see Eq. 12 
in [21]). For example, for VIq <C 1, the dominant met- 
ric mode decays as o^'^^q/^ during matter domination, 
which can also be obtained from Eqs. (|32ll33p . This in- 
troduces a significant suppression of large scale structure 
power during the matter era: 



5{nQ,z^Q) 
5{Q,z = 0) 



-30q/5 ^ o.6l"«/"-^ 



(37) 



if we fix the amplitude at matter-radiation equality, as- 
suming that it happened at z^q — 3500. Of course, only 
modes inside the horizon at matter-radiation equality will 
feel this amount of suppression. Modes entering after 
equality will suffer less suppression the later they enter. 
This leads to a red tilt in the cold dark matter spectrum 
up to the scale of matter radiation equality keq mimick- 



However, all scales 



ing a running spectral index [2 
k > keq are suppressed by the same factor ([37| which 
distinguishes quadratic Cuscuton from a simple running 
spectral index ACDM model. As the metric potential 
decays inside the horizon through the matter era, the 
CMB anisotropics receive an extra contribution from the 
Integrated Sachs- Wolfe (ISW). This leads to additional 
power, typically shifting the first Doppler peak of the 
CMB power spectrum towards slightly lower multipoles. 

Fig. ([2]) shows the resulting matter and CMB power 
spectra for VIq = 0,0.05, and 0.1, where we have fixed 
the initial amplitude of scale-invariant adiabatic pertur- 
bations, as well as the background expansion history, to 
that of WMAP3 concordance model Q. As expected, 
the changes in the amplitude of the matter power spec- 
trum (left panel in Fig. [5]) is most significant. However, 
direct measurement of the amplitude of the power spec- 
trum (e.g. from galaxy surveys) may be problematic due 
to the ambiguity in the value of linear bias. Methods to 
measure the linear bias, while present, are not as reliable 
as the shape of the power spectrum, as they involve non- 
linear physics, and thus are not widely used to obtain 
cosmological constraints. In addition, analytic calcula- 
tions of non-linear structure formation in similar model 
of early dark energy yield considerable and quite unex- 
pected deviations from the ACDM scenario precluding 
the use of standard approximations to infer the non-linear 
spectrum given the linear evolution [23j . 
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FIG. 2: Left: Dark matter power spectra for quadratic Cuscuton densities Qq — 0,0.05, 0.1, shown as solid, dotted, and dashed 
lines respectively. The initial amplitude of the scale-invariant adiabatic perturbations, is kept constant. Right: The CMB 
anisotropy power spectrum for the same cosmologies. 



Thus, instead of galaxy catalogs, we use Lyman-a for- 
est observations of quasar spectra from the Sloan Digital 
Sky Survey (SDSS) 5], which mainly constrains the lin- 
ear overdensity and spectral index ries at a pivot scale 
of fc ~ 0.009 s/km (~ 1 Mpc"^), and a pivot redshift of 
z ~ 3 In this mildly non-linear regime, we anticipate 
non-linear effects of to play little role. This point is 
further strengthened from the low values of Qq allowed 
by all other data without using Lyman-a (see below). At 
such low values of Qq, the remaining non-linear effects 
of fig will be negligible. Hence, Lyman-a forest observa- 
tions seem well suited to constrain the effects of fig on 
growth of linear perturbations. 

The changes in the CMB power spectrum are more 
subtle. As mentioned above, the main impacts of the 
Cuscuton quadratic potential can be seen in an addi- 
tional ISW contribution and a slight suppression of small 
scale power in the CMB spectrum. We ran a Monte 
Carlo Markov Chain (MCMC) simulation to constrain 
this model using a modified version of cmbeasv [23| . For 
this, we used the 3-year CMB data of WMAP "g. Super- 
novae la measurements constraints from b^ry- 
onic acoustic oscillations [27t. as well as the latest release 
of the SDSS galaxy survey [2^ and SDSS Lyman-a forest 
observations 5]. 

The WMAP 3-year data alone constrains fig < 3.6% 
at 95% confidence level. Adding SNe la and large scale 
structure data improves this - already substantial - 
bound further to Qg < 2.7%. Finally, adding Lyman- 
a forest data yields a tight limit of Qg < 1.6%. 

Hence, the background degeneracy of Ni, and ilg while 



present has little effect on the contours in Fig. The 
allowed abundance Qg is very small such that the effect 
on is at most AN^ = -3.04 x 2.47 x Qg x -0.2. 
This change in N^, is rather small compared to the 68% 
confidence interval for 7V^, which is Ni, = 5.5l5'2. 



C. Exponential (DGP-like) Potential 

DGP model of modified gravity Q, similar to other 
brane-world models, is inspired by non-perturbative ob- 
jects in string theory. It posits that the observable uni- 
verse lives on a 3-l-lD brane in 4+lD space-time. While 
matter fields are constrained to the brane, gravity can 
propagate into the bulk and thus can be sensitive to 
the full geometry of the Universe. The characteristic 
of the DGP model is that the brane action includes a 
term proportional to the volume integral of the induced 
Ricci scalar on our 3-brane. This is added to the or- 
dinary Einstein-Hilbert action in the bulk, which inte- 
grates the full Ricci scalar. It turns out that at high 
energies, the induced action dominates the dynamics of 
the induced metric, leading to ordinary four dimensional 
gravity, while at low energies, the bulk action takes over 
and gravity becomes five dimensional. This may lead to 
a self- accelerating phase of cosmic evolution, even in the 
absence of a cosmological constant. 

Although it is argued that the self-accelerating phase 
of the DGP model contains ghosts 31, 32, 33], and thus 
may not be realized in a physically stable way (but see 
[34|). it still remains the most widely studied concrete ex- 
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FIG. 3: Confidence contours for the abundance of the 
quadratic Cuscuton, and the number of relativistic neu- 
trino species A^,^ derived from using SNe la, CMB, large scale 
structure and Lyman-a forest observations. The white and 
black boundaries indicate 68% and 95% likelihood regions. 
While in principle, a degeneracy of the background evolution 
between SIq and Ni, exists, it is too weak to be seen here, be- 
cause the allowed abundance of Q.q is too small to have any 
significant effect on constraints (see the text). 



ample of a modified gravity model which competes with 
dark energy models as an explanation for the observed 
acceleration of the Universe. 

As we mentioned in Section [ill for a spatially flat 3- 
brane, the background evolution within self-accelerating 
DGP model coincides with that of a Cuscuton moAel with 
an exponential potential (Eq. [T2|) . Therefore, it is not 
possible to distinguish the two models based on the tra- 
ditional geometrical tests of background cosmology, such 
as supernovae la 0, Q, distance to the last scattering 
surface, or the scale of baryonic acoustic oscillations [13] ■ 

For sub- horizon perturbations, [sB'l and [29*1 compare 
the linear growth of structure in DGP model with a 
quintessence model with an identical expansion history. 
Since dark energy does not cluster on sub-horizon scales 
for Cs > 1, such a quintessence model should have a 
very similar behavior to the exponential Cuscuton model. 
[29| finds that the DGP model shows an extra 5-10% de- 
cay compared to the quintessence model (see their Fig. 
2). Therefore, it is not easy to distinguish DGP and 
a DGP-hke Cuscuton (or quintessence) model based on 
the growth of the large scale structure, [s^ estimate that 
this distinction can only be done at less 3cr level, with the 
next generation of supernovae, weak lensing, and (small 
angle) CMB observations. 

More significant is the ISW effect induced by the 
anisotropic stress in the DGP model, which is a char- 
acteristic feature of modified gravity models. [1^ find 
that, as a result of the ISW effect, the CMB anisotropy 
power spectrum is a factor of ~ 4 larger than than ACDM 
cosmology on small .^'s [1^. Therefore, they advocate 
use of correlations between CMB and high redshift large 




FIG. 4: CMB anisotropy power spectrum at low I's: The solid 
curves show the CMB power spectrum for the DGP-like dark 
energy models with = 1, 10, 100, where larger sound speeds 
have slightly larger power. Exponential Cuscuton is realized 
in the ct. ^ oo limit, although the power spectra for > 10 
are virtually indistinguishable. The dotted curve is the power 
spectrum for the ACDM cosmology, while the dashed region 
shows the prediction for the DGP self-accelerating model [291 ]. 
All the power spectra assume Q.m = 0.26 and h — 0.66. The 
points show the WMAP Syr observations, while the errorbars 
indicate an estimate of the cosmic variance error based on the 
data [M|. 



scale structure surveys (see e.g. i'^ order to de- 

tect signatures of the DGP model. As we show in Fig. 
([H), the CMB power spectrum for the exponential Cus- 
cuton model is nearly identical to that of the DGP-like 
quintessence model, and is still much smaller the predic- 
tion for the actual DGP model at low ^'s (47| . 

We thus conclude this section by pointing out that 
observational distinction between dark energy and mod- 
ified gravity models may be significantly more difficult 
than advocated in the literature (e.g. [Ill where the 
authors also point to a number of issues that need to 
be explored in such endeavor). A relatively simple and 
well-behaved dark energy model such as exponential Cus- 
cuton can exactly replicate the expansion history of the 
DGP self-accelerating cosmology, while predicting simi- 
lar small angle CMB and matter power spectra, at a few 
percent level. The smoking gun for the modified gravity 
models, thus, is their anisotropic stress, that can be po- 
tentially probed by the ISW effect in the low i regime of 
the CMB power spectrum. 
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V. ON THE CAUSALITY OF CUSCUTON 
LINEAR PERTURBATIONS 

In this section, we briefly comment on causal properties 
of the Cuscuton field theory with a minimal gravitational 
coupling. 

Eq. (|19p . which simply has a Yukawa screening form, 
can be written in the real space as: 



(5(^(x, t) — 3tpa^ 



-a|3ijli/^|x-x'| 

47r|x — x' 



-[$(x',i)+ifi>(x',t)]. 

(38) 

Therefore, we see that the non-local dependence of Cus- 
cuton on the metric perturbations is exponentially sup- 
pressed beyond a horizon that is defined by It is 
interesting to notice that an H horizon, also naturally 
occurs in cosmological gravitomagnetism, beyond which, 
the coupling of a gyroscope precession to the rotation of 
its surrounding environment saturates [s^ . 

Since metric scalar perturbations do not propagate, 
and tensor perturbations that do propagate, do not cou- 
ple to Cuscuton perturbations at the linear order, it is 
not possible to rigorously address the question of causal- 
ity of Cuscuton coupled to gravity, in the context of what 
we have done thus far. However, it is clear that the naive 
picture of instantaneous interaction, due to the infinite 
sound speed of Cuscuton, does not hold in a general rela- 
tivistic context, as Cuscuton is exponentially insensitive 
to metric fluctuations at separations larger than an H 
horizon. By the same token, and direct substitution into 
Eq. (|24p . one can see that super-horizon perturbations 
in the Bardeen parameter [40| : 

C = $ - + 

H 



= $ 



2(i/-i$ + $) 
3(1 + w) ' 



(39) 



remain constant until the modes enter the H horizon, 
where w is the effective equation of state. Therefore, 
there is no gross violation of causality in Cuscuton cos- 
mology. 



VI. CONCLUSIONS 

In this paper, we have studied the physical features, as 
well as possible observational signatures of a cosmology 
with Cuscuton dark energy, which was first introduced in 
a companion paper ^Ij, and could be realized as an in- 
compressible k-essence fluid. We showed that Cuscuton 
perturbations have no independent dynamical degree of 
freedom and, in lieu of other couplings, simply follow the 
space-time metric. Therefore, Cuscuton can be consid- 
ered to be a minimal theory of evolving dark energy, or 
a minimal modification of a cosmological constant. Due 
to lack of internal dynamics, Cuscuton only modifies (or 



dresses) the gravity of massive objects, and thus resem- 
bles a modified gravity theory. Indeed, to the best of 
our knowledge, this is the only modification of Einstein 
gravity that does not introduce any additional degree of 
freedom (and is not conformally equivalent to Einstein 
gravity). We then studied two specific Cuscuton cos- 
mologies, with quadratic and exponential potentials. 

We saw that the expansion history of Cuscuton cos- 
mology with a quadratic potential is identical to that of 
ACDM, and thus geometrical tests such as supernovae 
la, or the angular scale of barynoic acoustic oscillations, 
are blind to a quadratic term in the Cuscuton poten- 
tial. Nevertheless, the quadratic term acts as an early 
dark energy component with a constant energy fraction, 
and thus can be detected via its influence on the mat- 
ter power spectrum, or through the ISW effect in the 
CMB. We find that joint constraints from supernovae 
la, CMB anisotropics, power spectra of galaxy surveys, 
and Lyman-a forest fluctuations limit this component to 
fig < 1.6%, which is the most stringent upper limit that 
has ever been put on an early dark energy component. 

For the exponential Cuscuton model, we found that, 
surprisingly, the expansion history is identical to that of 
a flat DGP self-accelerating modified gravity model. The 
only observational distinction between the two cosmolo- 
gies is in the ISW effect at ^ < 20 in the CMB power 
spectrum, and is due to the anisotropic stress, present in 
the DGP modified gravity model. 

Indeed, cosmologists have yet to develop efficient tech- 
niques to detect a smoking gun for modified gravity. Ex- 
ponential Cuscuton cosmology is a clear example of where 
such a smoking gun may become necessary, and detec- 
tion of a non- vanishing anisotropic stress (see e.g. [HI), 
through ISW effect, appears to be the only way to rule 
out any such (scalar) dark energy model. An indepen- 
dent motivation for why new physics might be at work in 
the generation of large angle ISW effect may also come 
from observations of large angle anomalies in the CMB 
sky (e.g. [i^, and references therein). 

Finally, we showed that, despite its infinite speed of 
sound, Cuscuton perturbations are exponentially insen- 
sitive to the metric perturbations beyond the Hubble ra- 
dius, justifying why there is no gross violation of causality 
for superhorizon perturbations. 

It is now becoming clear that the phenomenology of 
dark energy will be the central theme in theoretical and 
observational cosmology, over the next ten years [43| . 
Nevertheless, all the current theoretical models for any 
observable deviation from a cosmological constant are 
at the best ill- motivated, or at the worst, already ruled 
out. Although falling short of solving the celebrated "cos- 
mological constant problem" f43| , Cuscuton dark energy 
provides an alternative which could be considered more 
natural, as it is protected against quantum corrections [l| 
(even in the presence of other couplings), and yet allows 
for a potentially observable evolving dark energy. The 
enthusiastic cosmologist may thus find some use for Cus- 
cuton in trying to convince the skeptical theorist about 
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the benefits of cosmological dark energy probes! 
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